ELECTRON-PHONON INTERACTION IN ULTRASMALL-RADIUS CARBON 

NANOTUBES 
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We perform analysis of the band structure, phonon dispersion, and electron-phonon interactions 
in three types of small-radius carbon nanotubes. We find that the (5,5) can be described well 
by the zone-folding method and the electron-phonon interaction is too small to support either a 
charge-density wave or superconductivity at realistic temperatures. For ultra-small (5,0) and (6,0) 
nanotubes we find that the large curvature makes these tubes metallic with a large density of states 
at the Fermi energy and leads to unusual electron-phonon interactions, with the dominant coupling 
coming from the out-of-plane phonon modes. By combining the frozen-phonon approximation with 
the RPA analysis of the giant Kohn anomaly in Id we find parameters of the effective Frolich 
Hamiltonian for the conduction electrons. Neglecting Coulomb interactions, we find that the (5,5) 
CNT remains stable to instabilities of the Fermi surface down to very low temperatures while for the 
(5,0) and (6,0) CNTs a CDW instability will occur. When we include a realistic model of Coulomb 
interaction we find that the charge-density wave remains dominant in the (6,0) CNT with Tcdw 
around 5 K while the charge- density wave instability is suppressed to very low temperatures in the 
(5,0) CNT, making superconductivity dominant with transition temperature around one Kelvin. 



I. INTRODUCTION 



It has been over a decade since the discovery of car- 
bon nanotubes (CNTs)i and the interest level in these 
systems continues to be high. The majority of theo- 
retical work on CNTs focuses on understanding the ef- 
fects of the electron-electron interactions using the cele- 
brated Luttinger liquid theory. 2 Experimental observa- 
tion of superconductivity in ropes of nanotubes^ and 
small-radius nanotubes in a zeolite matrix has also mo- 
tivated theoretical studies of the electron-phonon inter- 
actions (EPI), including the analysis of charge density 
wave (CDW)£i££& and superconducting (SC)V i n i 12 i 13 
instabilities. In this work we study the electron-phonon 
interactions in CNTs and discuss possible instabilities to 
the CDW and SC orders. Our approach provides reliable 
parameters for the effective Hamiltonians we use in con- 
trast to the Luttinger liquid treatments where obtaining 
such accurate quantities is quite difficult. 

A conventional starting point for discussing the 
electron-phonon interaction in solids is the Frolich 
Hamiltoniani^ 
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Here ct creates an electron with quasimomentum k in 
band r with spin er, a qfl creates a phonon with lattice 
momentum q and polarization /i, and q = k — k' mod- 
ulo a reciprocal lattice vector. The energies of electron 
quasiparticles and phonons (in the absence of EPC) are 
given by e kr and Q Q respectively. The EPC vertex is 
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Here \ip kr ) = c^ T |0) is a quasistationary electron state 
in band r with quasimomentum k, e qfi (i) is the phonon 
polarization vector on atom i in the unit cell, N c is the 
number of atoms per unit cell, M is the mass of a single C 
atom, N is the total number of unit cells in the system, 
and -jj^7 is the derivative of the crystal potential with 
respect to the ion position Roj. 

A common approach to obtaining parameters of the 
Hamiltonian Eq. for the CNTs is the zone-folding 
method (ZFM)»i£ The essence of this method is to take 
the electron band structure and the phonon dispersion 
for graphene and quantize momenta in the direction of 
the wrapping. The main results of such a procedure may 
be summarized as follows. The only bands crossing the 
Fermi level in graphene are the bonding and the anti- 
bonding combinations of the atomic p z orbitals. Hence, 
the zone-folding method predicts that these are the only 
bands which may cross the Fermi level in carbon nan- 
otubes. The condition for the quantized momenta to 
cross the Dirac points of the graphene gives the condi- 
tion for the (N,M) CNT to be metallic: N — M should be 
divisible by 3. The ZFM also predicts that the electron- 
phonon coupling in the CNTs should be dominated by 
the in-plane optical modes. This follows from the fact 
that the latter have the largest effect on the overlaps be- 
tween the p z orbitals of the neighboring carbon atoms. 

While the ZFM was shown to provide a quantitatively 
accurate description of the larger radius nanotubes, it 
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FIG. 1: The first Brillouin zone (a), electronic band structure 
(b), and phonon dispersion (c) of graphene. 



is expected to fail as the radius of the nanotubes is de- 
creased and the curvature of the C-C bonds becomes im- 
portant. Determining the band structure, the phonon 
dispersion, and the electron-phonon coupling of the small 
radius CNTs requires detailed microscopic calculations. 
In this paper we use the empirical tight-binding modeU^ 
to provide such analysis for three types of small-radius 
nanotubes: (5,0) with the diameter 3.9 A, (6,0) with the 
diameter 4.7 A, and (5,5) with the diameter 6.8 A. We 
find that the large curvature of the C-C bonds leads to 
qualitative changes in the band structure of the (5,0) and 
(6,0) nanotubes (band structures of the (6,0) and (7,0) 
nanotubes have been discussed previously ini£). For ex- 
ample, the (5,0) CNT becomes metallic from strong hy- 
bridization between the a and it bands (see Fig.^J). Fre- 
quencies of the phonon modes in small radius CNTs are 
also strongly renormalized from their values in graphene. 
Not only does the out-of-plane acoustic mode become a 
finite frequency breathing mode^ but even the optical 
modes change their energy appreciably (see e.g. Fig. CJ. 
Finally, the electron-phonon coupling changes qualita- 
tively in the small-radius CNTs. It is no longer dom- 
inated by the in-plane optical modes but by the out- 
of-plane optical modes which oscillate between the sp2 
bonding of graphene and the spy bonding of diamond (see 
discussion in Sec. IVIf) . We find that the strong effects of 
the CNT curvature decrease rapidly with increasing the 
tube radius. Already for the (5,5) nanotubes the ZFM 
gives a fairly accurate description of the band structure 
as well as the electron-phonon interactions. 

Determining parameters of the Frolich Hamiltonian for 
a one-dimensional system is not as straightfoward as for 
two and three-dimensional metals. Traditional methods 
for analyzing EPI from first-principles calculations are 
mean-field and, therefore, suffer from instabilities intrin- 



sic to one-dimensional systems. In particular, the frozen- 
phonon approximation, which is commonly used to de- 
termine the phonon frequencies, f^, in Eq. gives 
imaginary frequencies close to the nesting wave vector 
q = 2k f ■ This is the result of the giant Kohn anomaly, 18 
which corresponds to the Peierls instability of the one- 
dimensional electron-phonon system^ An important re- 
sult of our paper is that we developed a new formalism, 
which combines the frozen-phonon approximation with 
the Random-Phase Approximation (RPA) analysis of the 
EPI. This allows us to extract effective non-singular pa- 
rameters of the Frolich Hamiltonian from first-principles 
calculations or from the empirical tight-binding model. 
This technique should be applicable to many systems 
other than carbon nanotubes. 

After determining parameters of the Frolich Hamilto- 
nian Eq. for the (5,0), (6,0), and (5,5) CNTs we dis- 
cuss possible superconducting and charge-density wave 
instabilities in these systems. We find that neglecting 
the residual Coulomb interaction leads to much stronger 
CDW instabilities in all three cases (in such analysis 
Coulomb interaction is included only at the mean-field 
level via the energy of the single-particle quasi-stationary 
states, £fc r ). In the mean-field approximation we find 
the onset of the Peierls instability at temperatures 160, 
5, and 10~ 14 K for (5,0), (6,0), and (5,5) CNTs respec- 
tively. However, including the Coulomb interactions at 
the RPA level^S can lead to a stronger suppression of 
the CDW transition temperatures, TcdWj than the su- 
perconducting Tsc- For instance, we find by using the 
model Coulomb interaction of Ref. that for the (5,0) 
CNT, the CDW transition is suppressed to very low tem- 
peratures while superconductivity becomes the dominant 
phase with transition temperature of Tgc ~ 1 K. 

This paper is organized as follows. In Sec. [H] we dis- 
cuss our method for extracting parameters for the one- 
dimensional Frolich Hamiltonian. We then apply this 
method to the (5,0), (6,0), and (5,5) CNTs in Sec. UTTI 
In Sec. lIVI we use the constructed Hamiltonian for these 
CNTs to study their instabilities toward superconductiv- 
ity and charge-density wave states. The effect of intro- 
ducing the residual Coulomb interacting between elec- 
trons is covered in Sec. Finally all of the results are 
discussed and summarized in Sec. IV1I 



II. EXTRACTING PARAMETERS OF THE 
EFFECTIVE FROLICH HAMILTONIAN FROM 
THE FIRST PRINCIPLE CALCULATIONS 

Now we discuss our methods for calculating input pa- 
rameters to the Frolich Hamiltonian Eq. (JIJ for the rep- 
resentative nanotubes. Our analysis relies on the the em- 
pirical tight-binding modeU^ but it is easily amenable to 
any density-functional theory22i2^ treatment of the sys- 
tem. 
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A. Band Structure 

To compute the electronic structure of the CNTs we 
study, we use the NRL tight-binding method^ which has 
been tested and provided accurate results on a variety 
of materials. In this method, the Slater-Koster tight- 
binding matrix elements are parametrized and are fit 
to reproduce the first-principles density-functional band 
structures and total energies, with around 70 adjustable 
parameters per element. 

We study the (5,0), (6,0), and (5,5) CNTs which are 
shown in Figs. 1 1181 and The smallest possible unit 
cells for these CNTs contain 20, 24, and 20 atoms respec- 
tively. These CNTs are relaxed by minimizing their total 
energy per unit cell with respect to the atomic coordi- 
nates using 35 k-points in the first Brillouin zone. Matrix 
elements between neighboring atoms of up to 5.5 A were 
used. The calculations were performed on an orthorhom- 
bic lattice with spacing between parallel CNTs of 16 A, a 
distance sufficiently large to ensure negligible dispersion 
from inter-tube hopping. Once the CNTs are relaxed, 
the band structure is calculated. 



B. The phonon modes 

To calculate the electron-phonon coupling vertices and 
the phonon frequencies which will be discussed in the 
subsequent sections, one needs to have the ionic displace- 
ments corresponding to the normal vibrational modes of 
the CNT. As pointed out previously^ 7 we find that it 
is typically sufficient to use the zone-folded modes of a 
graphene sheet, even for the small-radius CNTs we study 
as will be discussed below. 

Following the method used in the book of Saito et al.^ 
we have computed the 60 x 60 dynamical matrix of a (5,0) 
CNT and in Fig. [21 we compare the resulting phonon dis- 
persions with the zone- folding results. The ionic displace- 
ment modes obtained by the two different methods are 
very similar except for a few special cases. For instance, 
the zone-folding results give three acoustic modes which 
correspond to translating the graphene sheet in differ- 
ent directions. Upon rolling the graphene sheet, these 
modes get mapped to two acoustic modes corresponding 
to rotation about the CNT axis and translation along 
the CNT axis and the optical breathing mode. Con- 
versely, diagonalizing the dynamical matrix of the CNT 
gives four acoustic modes corresponding to translations 
in three directions and the rotating mode (actually us- 
ing the method of Ref. Il5l one obtains a small spurious 
frequency for the rotating mode as pointed in this refer- 
ence). Upon unrolling the CNT to the graphene sheet, 
the rotating mode and the mode corresponding to trans- 
lation along the CNT axis will become acoustic modes of 
the graphene sheet. However, the two CNT translational 
modes which are perpendicular to the CNT axis will get 
mapped to ionic displacements which are not eigenmodes 
of the graphene sheet which are mixtures of in-plane and 
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FIG. 2: The phonon dispersions of a (5,0) CNT determined 
by (a) the zone-folding method and (b) diagonalizing the full 
dynamical matrix of the CNT. 

out-of-plane oscillations. In addition, using the dynami- 
cal matrix of the CNT, we find that there is mixing be- 
tween the breathing and strething modes around k = 0.3. 
In this vicinity, there is level repulsion from the lifting of 
the degeneracy of these modes. Away from this point, 
the modes are, to a good approximation, decoupled. 

In our analysis of the electron-phonon coupling we use 
the displacements obtained from the zone-folding method 
to simplify the calculations, as well as to give a clear 
conceptual picture. We then check that none of the im- 
portant electron-phonon couplings come from any of the 
few graphene modes for which the zone-folding method 
breaks down. 



C. The electron-phonon coupling vertices 

The electron-phonon coupling (EPC) matrix in Eq. © 
can be evaluated by using the finite difference formula 

AffaW/. = \{>Pkr\{V qli ~ V )\^k'r') (4) 

where V qfl and Vq are the perturbed and the unperturbed 
lattice potentials respectively. A method for calculating 
the expression with a plane-wave basis set was pre- 
viously developed^* In this paper we extend this proce- 
dure to tight-binding models. We introduce the standard 
tight-binding notation 

\lpkr) = ^2 A kTl l\xkil) ■ (5) 

ll 



\xku) = -±=J2 eik ' Rn ^- ( 6 ) 

* n 



Here \<f> n i) are the electron states for isolated carbon 
atoms, n runs over unit cells, i runs over basis vectors 
in the unit cell, and I runs over orbital type. We find (for 
details, see Appendix 

A/ fcTfe , T> = I £ AUdxZm™ - E F )\x q k ^ v )A klrHn , 

ili'V 

(7) 

This expression can be computed by evaluating the tight- 
binding Hamiltonian and overlap matrices for the dis- 
torted lattice, evaluating the coefficients and Ayy of 
the wave functions for the undistorted lattice, and per- 
forming the above sum. 

In all the calculations presented in this paper we used 
the ZFM to find phonon eigenvectors in the nanotubes 
starting from the phonon eigenvectors in graphene. 15 The 
latter have been obtained using the 6x6 dynamical ma- 
trix of graphene given in Ref. |25[ We emphasize that 
we use the ZFM only to find the phonon eigenvectors in 
small nanotubes, but not the phonon frequencies. The 
frequencies are affected strongly by the CNT curvature, 
and should be computed directly. This is discussed in 
detail in Sec. HTTJ1 and Sec.lITTI 



D. Phonon frequencies 

A standard method of calculating the bare phonon fre- 
quencies Cl® in Eq. is the frozen-phonon approxima- 
tion (FPA)i2& In this approach 



i 



:y/AE cos (q) + AE sin (q) 



(8) 



where u is the amplitude of the displacement, and 
AE cos (q) and AE s [ n (q) are the energy differences per 
unit cell between the distorted and equilibrium lattice 
structures where the distortion corresponds to the real 
and imaginary parts of SH n i = ue tqR ™e qiJ _(i) respec- 
tively. When we apply this procedure to one-dimensional 
CNTs, we find that AE cos (q) + AE S i n (q) becomes neg- 
ative around certain wave vectors (see e.g. Fig. [7J. A 
closer inspection shows that such anomalous softening 
always corresponds to one of the 2k F wave vectors of the 
electron bands indicating the presence of the giant Kohn 
anomaly. 

It is important to realize that the divergence of Q qfl 
obtained in the FPA does not imply the divergence of 
in the Frolich Hamiltonian Eq. (JJJ. The frequencies 
Q qfJ , are calculated after the electron-phonon interaction 
in Eq. QJ have been included, which gives anomalous 
softening at 2kp due to the well-known Peierls instabil- 
ity of electron-phonon systems in Id. In two and three 
dimensional systems renormalization of the phonon fre- 
quency by electrons in the conduction band is typically 
negligible. So, one can use phonon energies obtained in 
the FPA as a direct input into the Frolich Hamiltonian. 
By contrast, nesting of the one-dimensional Fermi sur- 



FIG. 3: The phonon propagator evaluated within the RPA. 



faces, leads to dramatic renormalization of the phonon 
dispersion by electrons in the conduction band. 

To extract the bare phonon frequency f2° from the 
numerically computed we point out a connection 

between the FPA and the RPA for the Frolich Hamilto- 
nian. For negligible interband coupling (this condition is 
satisfied for all modes showing the giant Kohn anomaly, 
which we discuss in this paper) Dyson's equation for the 
phonon propagator _D(g,ii/ m ), as shown in [21 is given by 

D^iUm) = D 0fl (q,iu m ) (1 +Jl, J ,(q,iv m )D lt (q,iu m )) . 

(9) 

Here v m = 2nmT are the bosonic Matsubara frequencies 
and 



D 0fl (q,iu m ) = 



2n 



(10) 



is the non-interacting phonon Green's function. The 
phonon self-energy evaluated in the RPA is given by 

IL^(q,ii/ m ) = 2T V" \g pTp+qT ^\ 2 Gor(p+q, iio m+n )G 0T (p, iw n ). 



npr 



where non-interacting electronic Green's functions are 
given by Gq t (p, iu> n ) = (iujn—epr)^ 1 andw n = 7r(2n+l)T 
for integer n are the fermionic Matsubara frequencies. 
Summing over n, we obtain for Eq. (|ll|l 



where the bare susceptibility is given by 

f{s pT ) — fi^p+qr) 



£pr £p-\-qr 



(12) 



(13) 



with fispr) = (1 + e l3e i> T )~ 1 being the Fermi-Dirac dis- 
tribution function. 

The poles of the phonon Green's function £l qp (we put 
iv m — * Qqp, in D p (q, iv m )), which give the dressed phonon 
frequencies, will satisfy the equation 

WSpa = W 2 + ^^(q, ««.)• (14) 

Due to the large energy difference between electrons and 
phonons, it is typically a good approximation to set 



a, 



in n M (g, Qqfj,). This approximation results in 



an expression that can be derived by doing stationary 
second-order perturbation theory to obtain the change 
in energy due to the presence of the phonon. That is, 
setting £l ql± — > in H(q, Q qt i) corresponds to the frozen- 
phonon approximation 



(*V) 



WFPA 



■2n° n M ( 9) o). 



(15) 



5 




FIG. 4: The band structure of the (5,0) CNT obtained 
through zone-folding (a) and calculated directly (b) along 
with the atomic structure (c). 



We can typically approximate well the quasiparticle en- 
ergy by a plane- wave state with given effective mass m* . 
Then, by incorporating the FPA, at zero temperature the 
integral in Eq. 111. 'ill can be done which will enable us to 
obtain 
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(16) 



This expression explicitly shows the logarithmic diver- 
gences in the phonon dispersion at the nesting wave vec- 
tors of the Fermi surface. This is the famous Peierls in- 
stability to a CDW state. Our procedure for determining 
the elusive undressed frequencies is then as follows. We 
take f2 9At obtained from the FPA and fit them with the 
expression Eq. lfTo|l using fi° as an adjustable parame- 
ter. The coefficients of the log divergences at the nesting 
wave vectors of the Fermi surface are fixed by the effec- 
tive masses m* and k Fr (known from the band structure) 
and the computed EPC matrix elements M 2 fc FT/J . In all 
cases we found excellent agreement of the calculated FPA 
frequencies with Eq. 1|16|) in the vicinity of the singular 
points, which provides a good self-consistency check for 
our analysis. 



III. RESULTS FOR REPRESENTATIVE 
NANOTUBES 

A. (5,0) nanotube 

The zig-zag (5,0) CNT has a diameter of around 3.9 
A making it close to the theoretical limit i^i Nanotubes 
of this size have been experimentally realized through 
growth in the channels of a zeolite host. 4 Through the Ra- 
man measurement of the frequency of the radial breath- 
ing mode, the (5,0) CNT is thought to be a likely candi- 
date structure for these experiments. 28 

We first compute the band structure of this tube 
by using the zone-folding method^ To do this, we 



use the band structure of graphene, which is shown in 
Fig-ffl computed by using the NRL tight binding method. 
Shown in this figure are four valence bands and four con- 
duction bands, coming from the three spi and one p z 
bonding and antibonding states respectively. There is 
a degeneracy between the p z bonding and antibonding 
states at the Fermi energy at the K point in the first 
Brillouin zone which accounts for the semimetallic be- 
havior of graphene. The zone-folding band structure of 
the (5,0) CNT is shown in the right of Fig. H Since 5/3 
is not an, zone-folding predicts this CNT to be semicon- 
ducting. 

Fig.H(b) shows the band structure of the (5,0) CNT 
calculated directly by using a unit cell of 20 atoms. One 
sees that there are significant qualitative differences be- 
tween the two band structures, one being that the di- 
rectly computed band structure predicts metallic behav- 
ior. The inner band (with smaller Fermi point kp) is dou- 
bly degenerate while the outer band (with larger Fermi 
point kp) is nondegenerate. The strong curvature effects 
causes hybridization between a and ir bands, pushing 
them through the Fermi energy and therefore making 
the tube metallic. Furthermore, for the (5,0) CNT, we 
see that inner band is close to the Van Hove singular- 
ity at k = 0, which produces a large density of states at 
the Fermi energy. The calculated density of ^(0) = 0.16 
states/ eV / C atom is around a factor of five larger than 
that of larger radius metallic armchair CNTs. 

After the band structure is calculated, we consider all 
possible scattering processes of electrons between Fermi 
points —kp,—kp,kp, and k F due to phonons with wave 
vectors q that satisfy the momentum conservation con- 
dition. As a starting point for the phonon spectrum, we 
use the dynamical matrix of Jishi et alm^ which uses a 
fourth nearest-neighbor model, and we employ the zone- 
folding method. The reproduced phonon dispersion of 
graphene is shown in Fig. ^ For a given process, we cal- 
culate the coupling for all of the 3 x N c distinct phonon 
modes where N c — 20 is the number of atoms per unit 
cell. Shown in Fig.JSJis an example of the outcome for one 
of these calculations. Shown is the coupling for the outer 
band 2kp processes vs. graphene frequency. One can im- 
mediately see that most couplings vanish which can be 
explained by symmetry of the electronic wave functions 
and the phonon modes. 

To keep this paper concise, we cannot present all of 
the coupling results for each scattering process. Instead, 
we show the most dominant couplings. These dominant 
couplings were found to be from intraband 2kp pro- 
cesses. The largest couplings for the (5,0) CNT occur 
for phonons along the TM line of graphene at the appro- 
priate wave vector corresponding to the particular 2k p. 
For the inner band, the largest couplings, in descend- 
ing order, occur for the out-of-plane optical mode, the 
radial breathing mode, and the in-plane acoustic stretch- 
ing mode. For the outer (with larger kp) band, the dom- 
inant couplings occur for the out-of-plane optical, an in- 
plane optical, the radial breathing, and in-plane stretch- 
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FIG. 5: The coupling M^y for the outer band 2k F process 
for each of the 3 x 20 = 60 phonon modes respectively vs. 
graphene frequency. 
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4.10 



TABLE I: Calculated values for the dominant coupling pro- 
cesses for the (5,0) CNT. The numbering scheme here corre- 
sponds to that given in Fig. |S| 2kp and 2k p correspond to 
inner and outer band processes respectively. Phonon frequen- 
cies are given for graphene. 




(a) 



FIG. 7: (a): Phonon dispersion for the (5,0) CNT along the 
TM line of graphene. The X's denote values for which the 
frozen-phonon approximation gave imaginary frequencies for 
the out-of-plane optical mode in the vicinity of 2kp. (b): The 
mode showing the most softening fit to the RPA expression. 



FIG. 6: 1-6: The phonon modes at the F point in the first 
Brillouin zone of graphene. 7: An in-plane optical phonon 
mode at the K point of the first Brillouin zone of graphene. 
The out-of-plane optical mode 4 is the leading cause of the 
CDW instability in the (5,0) and (6,0) CNTs. 



ing modes. These results are summarized in Fig. and 
Table H| Although the magnitunde of the dominant cou- 
pling matrix element for the outer band is larger than 
that of the inner band, the inner band processes are sig- 
nificantly more important in the study of instabilities be- 
cause their contribution to the total density of states at 
the Fermi energy is significantly larger than that of the 
outer band. This is due to the small Fermi velocity of 
the inner band and its degeneracy. 

It is interesting to note that the phonons that have the 
strongest coupling to electrons at the Fermi surface are 
out-of-plane modes. This is different than intercalated 
graphene where in-plane phonon modes are responsible 
for superconductivity^ The fact that the out-of-plane 
modes are the most important for this CNT are presum- 
ably due to the large curvature effects. For instance, we 
find that the bond angles of the relaxed (5,0) CNT struc- 
ture (having the values of 119.4° and 111.9°)are interme- 
diate between the sp2 bond angle (found in graphene) 
of 120° and the sp3 bond angles (found in diamond) of 
109.4°. 



Now we calculate the CNT phonon frequencies by us- 
ing the frozen-phonon approximation with the eigenvec- 
tors from graphene. The circles shown in Fig. \7\ are the 
frequencies obtained for phonon modes along the TM 
line of graphene for the out-of-plane optical mode which 
was found to be the most important mode. First, we 
see that the calculated FPA frequencies are significantly 
lower than the corresponding ones in graphene. This can 
be understood as follows. The strong curvature of the 
nanotube changes the C-C bonds so that they are in an 
intermediate regime between the sp2 bonding (found in 
graphene) and sp^ bonding (found in diamond) . The out- 
of-plane optical mode oscillates between these two bond- 
ing configurations and is therefore significantly softened. 
Next, we notice that there are divergences at q = 2k p 
and q = 2kp. This result is the giant Kohn anomaly. 

To extract the bare phonon frequency of the Frolich 
Hamiltonian Eq. 2] for the (5,0) CNT we follow the pro- 
cedure discussed in Sec. Ill Dl The dressed phonon fre- 
quencies are given by 
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All of the quantities needed to calculate the coefficients 
Z?A and Dehave been obtained already. We assume that 
the bare phonon frequencies are fit well by the form 
(f2° M ) 2 = ao + a\q + a2q 2 ■ We then use <2o,<Xi and a<x 
as fitting parameters to fit our expression for f2 g to the 
calculated FPA frequencies. Doing this thereby enables 
us to extract the important bare frequency dispersion 
which is shown in Fig. [7\ Extracting these bare fre- 
quencies Sl^ allows us to calibrate the effective Frolich 
Hamiltonian Eq. Q which will be used to study insta- 
bilities of the electron-phonon system. With our previ- 
ously calculated quantities, we obtain Da = (219 cm -1 ) 2 
and Db = (146 cm" 1 ) 2 . Using these values we thereby 
extract il° „, A — 433 cm -1 . 



B. (6,0) nanotube 

The band structure of the (6,0) CNT was considered 
extensively by Blase et al. in Ref. 0. This tube has a 
slightly larger diameter of 4.7 A. The zone-folding band 
structure of this CNT is shown in the left of Fig. [S] As 
is typical of metallic zig-zag tubes, there are two bands 
crossing at k — at the Fermi energy. The band struc- 
ture directly computed with 24 atoms in the unit cell 
is shown in the right of Fig. [S] As discussed beforeii, 
these band structures differ qualitatively which is a re- 
sult of the hybridization of the sp2 and p z bands. Here 
the inner band (with smaller k F ) is nondegenerate and 
originates from the p z bonds in graphene while the outer 
band (with larger kp) is degenerate and originates from 
the sp2 bonds in graphene. 

The coupling matrix elements for the (6,0) CNT were 
computed and the coupling for the most dominant modes 
are shown in Fig. and Table [n] The dominant inner 
band couplings were for intraband processes and are, in 
descending order, to the out-of-plane optical and an in- 
plane optical. The dominant outer band couplings pro- 
cesses were found to be the out-of-plane optical mode, 
an in-plane optical mode, the radial breathing mode, and 
the in-plane acoustic stretching mode. 

Using the same procedure as was used for the (5,0) 
CNT in the previous section for extracting the bare 
phonon frequency at 2k p. From the previously computed 
values for the electron-phonon coupling matrix elements 
and the band structure, we find Da = (166 cm^ 1 ) 2 and 
Db = (107 cm -1 ) 2 . After fitting, we extract the value 
^£=480 cm" 1 . 




FIG. 8: The band structure of the (6,0) CNT obtained 
through zone-folding (a) and calculated directly (b) along 
with the atomic structure (c). 



(6,0) 


mode 


wf aph (cm" 1 ) 


M kk ,(eV/A) 




4 


857 


7.27 




5 


1585 


6.80 




4 


847 


6.84 




6 


1591 


6.12 




3 


68 


3.73 




2 


493 


2.31 



TABLE II: Calculated values for the dominant coupling pro- 
cesses for the (6,0) CNT. The numbering scheme here corre- 
sponds to that given in Fig. |S| 2kp and 2fc^ correspond to 
inner and outer band processes respectively. Phonon frequen- 
cies are given for graphene. 



C. (5,5) nanotube 

Finally, we study the more conventional armchair (5,5) 
CNT which has a diameter of around 6.8 A. As shown 
in Fig. the zone-folding and directly computed band 
structure for this larger diameter tube agree quite will. 
Both of these band structures show two bands which orig- 
inate from p z orbitals which cross at the Fermi energy at 
around k = | - . 

3 7T 

The largest couplings for the CNT were found to again 
be from the intraband processes and are shown in Fig. 
and Table ITTll The only significant intraband coupling is 




(a) 



FIG. 9: (a): Phonon dispersion for the (6,0) CNT along the 
TM line of graphene. (b) : The mode showing the most soft- 
ening fit to the RPA expression. 
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FIG. 10: The band structure of the (5,5) CNT obtained 
through zone-folding and calculated directly are shown in the 
upper right and left. Bottom: the structure of the (5,5) CNT. 



(5,5) 


mode 


tu 3 q ravh (cm" 1 ) 


M fcfc ,(eV/A) 




7 


1479 


11.60 




3 


542 


4.64 



TABLE III: Calculated values for the dominant coupling pro- 
cesses for the (5,5) CNT. The numbering scheme here corre- 
sponds to that given in Fig. |S| Phonon frequencies are given 
for graphene. 



for an in-plane mode shown denoted by 7 in Fig. [S] The 
wave vector for this mode is at the K point in the first 
Brillouin zone of graphene. For the interband processes, 
there is coupling to the the radial breathing mode, but 
this is significantly smaller. 

For the (5,5) CNT, applying our method of extract- 
ing the bare phonon frequencies, we obtain Db = (228 
cm -1 ) 2 . Note that for this system, only it bands are rel- 
evant at the Fermi surface. We extract Sl° _, B = 1469 

cm -1 . 

It is worth pointing out that there has been some con- 
troversy about the relevant phonon mode which couples 
the electrons at the Fermi surface for the (5,5) CNT&2S 
Our results confirm the study of Ref. 13(1 The 2fcp pro- 
cesses couple to the phonons at the K point of graphene 
and the relevant graphene mode has polarization vectors 
e,(l) = 775(2,1, 0) and e,(2) = ^g(l,t,0). This out-of- 




(a) 



FIG. 11: (a): Phonon dispersion for the (5,5) CNT along 
the MK line of graphene. (b): The mode showing the most 
softening fit to the RPA expression. 





(5,0) 


(6,0) 


(5,5) 


mode 


4 


4 


7 




433 


480 


1469 


Acdw 


0.26 


0.12 


0.024 


Tcdw (K) 


160 


5 


7 x 10" 14 



TABLE IV: The dominant mode for the CDW instability, the 
extracted bare phonon frequency, the CDW coupling param- 
eter, and the CDW transition temperature for the various 
CNTs studied. 



phase circular motion is qualitatively different from the 
linear oscillations thought to couple previously. 



IV. INSTABILITIES OF THE 
ELECTRON-PHONON SYSTEM 

A. Charge-Density Wave Order 

The RPA analysis presented in Sec. Ill Dl can be used 
to investigate the CDW (Peierls) transition temperature. 
This instability corresponds to softening of the phonon 
frequency to zero, so we can obtain it from the condi- 
tion Qq t = in Eq. (|14fl where Q T = 2k-p T is one of 
the nesting wave vectors of the Fermi surface. The elec- 
tron polarization evaluated at temperature T is given 
by X 0r(2k FT ,uj = Q,T) = |i/ T (0) log(T/4e Fr ), where 
i/ T (0) = 2m* / 'Lk-pr is the contribution to the total den- 
sity of states from band r. We introduce the CDW cou- 
pling constant 



Ac 



DW,r/i 



\9Q^\ 2 M0) 

n° . 



(20) 



where r specifies which of the 2k p nesting wave vectors 
we are considering and fi labels the phonon mode. Note, 
that distinguishing between various phonon modes is im- 
portant, since it tells us about the nature of the distortion 
of atoms below the Peierls transition (i.e. the in the plane 
vs out of the plane) . One finds for the CDW transition 
temperature 



CDW,T/i 



= 4£ 



-1/Acdwt/j 



(21) 



Corrections to this equation due to an additional band 
with different Fermi wave vector (e.g. the term with the 
logarithmic divergence at 2k J? in Ea. 117(1 is small and will 
be neglected. Degenerate bands (e.g. the A band for the 
(5,0) CNT), are accounted for by an additional factor of 
2 in the density of states is Eq. J3UJ. In Table llVl W e 
summarize our results for the CDW instability for the 
CNTs studied. 



B. Superconductivity 

To analyze the superconducting instability of the 
CNTs we use the Migdal-Eliashberg theory. The 
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isotropic Eliashberg equations for the one-dimensional 
case, neglecting the Coulomb interaction, can be writ- 
ten as (see Appendix IbI for details) 



Z n 



1 + f n S n ^ K n - n ') s n' 



Z n Ar 



E 



A(n - n')/„/A r 



(22) 



(23) 



where /„ = l/\2n+ 1|, s„ = sgn(2n + 1), A„ = <j) n /Z n , 
and the frequency dependent coupling constant X(n) is 
given by 



X(n — n) 



1 



^(0) 



Eferfc'r> 5(£ k r)5{e klT ')\gkk'^\ 2 

X D^{k-k',n-n'){2A) 



where v a (0) is the density of states per spin at the Fermi 
energy. When analyzing superconductivity in two and 
three dimensional systems using the Eliashberg equa- 
tions it is sufficient to take the bare phonon propagators 
Do(fc — k' , n — n') in Eq. I|24|l . This is justified since in the 
absence of Fermi surface nesting there is typically little 
difference between the bare and the dressed phonon fre- 
quencies and propagators. In one-dimensional systems, 
however, there is a strong temperature dependent renor- 
malization of the phonon spectrum which needs to be 
taken into account. The simplest way to do so is to use 
the FPA form of the phonon propagator (see eqs © - 



2n° 



(iv m ) 2 - (^V) 2 



(25) 



Here f2 9AI is the dressed phonon frequency in the FPA 
given in Eq. (fT5|l . Taking a soft dressed phonon propa- 
gator immediately leads to the enhancement of the elec- 
tron pairing via the increase of A(n). Enhancement of 
superconductivity by the giant Kohn anomaly in onc- 
dimcnsional systems has been discussed previously by 
Heeger in Ref. 31. The main subtlety of the Eliash- 
berg equations in this case is that the phonon frequency 
Q ql/ now has temperature dependence which needs to be 
found using the finite temperature form of the polariza- 
tion operator U(q, 0) in Eq. (tT3|) . 

When we analyze the (5,0) nanotube following this 
strategy, we find, however, that the CDW instability al- 
ways appears before the superconducting one. This is in 
agreement with the general argument proposed in Ref. 32 
that in strictly one-dimensional electron-phonon systems 
Peierls instability alway dominates, since it involves all 
electrons in the band, compared to the superconducting 
instability, which involves only electrons in the vicinity 
of the Fermi surface. 

To introduce a quantitative measure of the strength of 
superconducting pairing we use the bare phonon prop- 
agator in Eq. I|24|) . This approximation will be more 





(5,0) 


(6,0) 


(5,5) 


Asc 


0.57 


0.12 


0.031 


r sc (K) 


64 


0.071 


1.11 x 10" 12 



TABLE V: The SC coupling parameter, and the SC transition 
temperature for the various CNTs studied. The CDW insta- 
bility and the residual Coulomb interaction between electrons 
are neglected in the calculation of these quantities. 



carefully considered in Sec. IV Dl along with inclusion 
of the Coulomb interaction. A useful approximate so- 
lution of the Eliashberg equations (1221) - l|24|l is given by 
the McMillan formula (again in the absence of Coulomb 
interaction)2ii2i 



Tsc 



1.20 



exp 



1.04(1 + A 



sc) 



Asc 



(26) 



Here Asc is the zero frequency component of Eq. I|24|) 
where, again, the bare phonon frequencies are used 



Asc 



1 



MO) 



J2 5{e k )5(e k ,)\g kkl \ 2 D (k - k' , 0). (27) 



kk' 



In accordance with Ref. H we take (fi) = 1400 K. The 
superconducting coupling constants and transition tem- 
peratures calculated in this manner are summarized in 
Table E] We emphasize, however, that these numbers 
should be taken with some scepticism, since within the 
same approximation the CDW instability is usually the 
dominant one and appears at much higher temperatures 
(compare to Table Hvjl . 

Finally, it is known that q sa scattering processes due 
to acoustic phonons can be important in one-dimensional 
electron-phonon systems £££§£1 However, in the approx- 
imations leading to Eq. 1271 these contributions were ne- 
glected. In Appendix El we show that while these pro- 
cesses can be important for some systems, their inclusion 
leads to only a small correction to Age for the CNTs we 
study. This is due to the fact that the dominant con- 
tributions to the superconducting coupling constant are 
from optical phonons. 



V. ROLE OF THE COULOMB INTERACTION 

In the discussion above we concentrated on the 
electron-phonon interaction with electron-electron 
Coulomb interaction included only at the mean-field 
level via the band structure. It is useful to consider 
how the residual Coulomb interaction can modify the 
analysis of the Peierls and superconducting instabilities 
discussed above. We take 



H — 7i -ph + W c -c (28) 

= 2 V<lTT'C k +qTcr C k'-qT'cr' C k'T'a<C k Tcr 



kk' qrr' aa' 
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where Tt c - P h is still given by Eq. (JTJ and we will always 
assume k and k' around the Fermi surface. Note that 
we have neglected interband scattering which is typically 
small. In the following, we will consider how introducing 
this Coulomb interaction modifies the results. 



A. Coulomb interaction potential 



(a) 



(b) 



FIG. 12: Dyson's equation for the phonon propagator (a) 
where the Coulomb interactions are taken into account within 
the RPA (b). 



For the Coulomb interaction between conduction elec- 
trons, we take the form used by Egger et al. in Ref. 



B. Modification of CDW instability due to 
Coulomb interaction 



V(r-r') = 



e 2 / K 



(x-x') 2 +[2Rsm[lijf)) + a 2 



(29) 

Here, the y-direction is chosen to be along the perimeter 
of the CNT and x measures the distance along the CNT 
axis. A measure of the spatial extent of the p z electrons 
perpendicular to the CNT is given by a z w 1.6 A and R 
is the CNT radius. Note that this interaction potential 
is periodic in the y-direction. For the dielectric constant 
due to the bound electrons, we will take the value k rs 2 
predicted by the model of Ref. 

We can now use Eq. (|29|) to obtain the Coulomb inter- 
action entering Eq. lj2l3|) 



The simplest approximation (beyond mean-field) 
which includes the Coulo mb rep ulsion is the RPA shown 
in Fig. EKsee e.g. Refs. Il4l20lh Eq. (jnj now becomes 
for a one-band system 



'1-V(g)xd(«,n, 



(32) 



where U M (q,Sl qfl ) = \g q ^\ 2 Xo(q, tt qfl ). We immediately 
see that including the Coulomb interaction can suppress 
the CDW instability. The second term in Eq. (|32[l no 
longer diverges when q = 2kp and the softening of the 
2kp phonons occurs only for hcdw < Aqdw ui where 



V qTT , = J d 2 rd 2 r'V(r - r') (30) 

The region of integration above is over areas of length L 
along the cc-direction where L is the length of the system 
and of width 2irR along the y-direction. For backward 
scattering processes (q » 2kp) between the inner bands 
of the (5,0) and (6,0) CNTs we find that V qTT ' is inde- 
pendent of r and t', and (see Appendix El for derivation) 



V q » 7-^2 / dxdx'e- tq{x ' x,) 



2ttR 



dy 
2nR 



2irR 



dy' 
2-kR 



(31) 



where 7 = 0.59 and 0.0016 for the (5,0) and (6,0) CNTs 
respectively. This is significantly reduced from the value 
of 7 w 1 that one obtains for larger radius CNTs^i which 
is due to the fact that wave functions at the Fermi points 
have different symmetries for the (5,0) CNT and (6,0) 
CNTs. More specifically, it can be found that /icdw is 
very small for the (6,0) CNT due to the fact that for 
metallic zig-zag nanotubes, the wave functions at —k 
and k close to the Fermi energy are nearly orthogonal 
within the unit cell of the CNT since they correspond 
to symmetric and antisymmetric combinations of atomic 
orbitals in the graphene sheet. 



Mcdw = ~v(0)V q =2ki 



(33) 



From equation l/i2[l we also find how the Coulomb inter- 
action modifies the Peierls transition temperature 



Tcdw.ju = 4£f exp 



[A 



CDW.p — MCDWj 



(34) 



We will now estimate the magnitude of /xcdw from 
this residual Coulomb interaction for the (5,0) which was 
seen above to be the most unstable toward the forma- 
tion of a CDW from distortion of the out-of-plane optical 
mode shown in Fig. EI Carrying through the straightfor- 
ward generalization of the RPA analysis for the multiple- 
band system, and carrying out the integrals in Eq. I|31|) 
for the Coulomb backward scattering interaction, we ob- 
tain hcdw = 0.24. Note that this is quite close to 
Acdw = 0-26 for this particular instability. This indi- 
cates that it is possible that the Coulomb interaction can 
significantly suppress the CDW transition temperature 
or even remove the CDW instability altogether. Indeed, 
taking these values we find that Tcdw is suppressed to 
less than 10" 18 K. 

For the (6,0) CNT, we calculate the smaller value 
Mcdw = 0.0019. This will not change the value of 
Tcdw = 5-0 K that we calculated previously for the (6,0) 
CNT. 
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Phonon vertex renormalization through 
screening 



find that the renormalized contribution to the supercon- 
ducting coupling constant is given by 



It can be seen that the Coulomb interaction further can 
screen the electron-phonon vertex. By including screen- 
ing through the RPA, we find that the screened vertex is 
given byi* 



9q» 



1 - V qX o(q) 



(35) 



for a one-band system where Xo(<?) = Xo(<7, ^<?m = *-*)• 
Thus we see that the inclusion of screening reduces the 
electron-phonon vertex. We note that in the treatment 
in Sec. IV Bl of the CDW instability it would be inappro- 
priate to use the screened vertices since this would lead 
to double-counting. 

Full charge self-consistent calculations will determine 
the dressed electron-phonon vertex (see Appendices E] 
and^J). This is desirable in 3d, where the renormaliza- 
tion is presumably small. However in Id, one would cal- 
culate greatly suppressed values for the couplings, domi- 
nated by the screening due to the logarithmic divergence 
of the susceptibility at 2kp . Because of the subtle inter- 
play between these divergences, it is desirable to calculate 
the bare vertex and then manually put in the Coulomb 
interaction as we do. 

Since with the method we use, the charge distribu- 
tion is not calculated self-consistently, we calculate the 
bare electron-phonon vertex g qfl . We point out, how- 
ever, that there is an approximation here. The true bare 
electron-phonon vertex should be calculated in the ab- 
sence of the conduction electron entirely which is sepa- 
rately accounted for in the residual Coulomb term. In our 
method, however, the conduction electron is taken to adi- 
abatically follow the ion through the distortion. Because 
of this, we expect our results to slightly underestimate 
the bare electron-phonon coupling vertices. 



D. Modification of superconducting instability due 
to Coulomb interactions 



To include the Coulomb interaction in the Eliashberg 
equations, it is necessary to dress both electron-phonon 
vertices shown in Fig. 1141 according to Sec. IV"Cl as well 
as the phonon propagator according to Sec. IV Bl This 
leads to the modified phonon-mediated interaction be- 
tween electrons of 



\9qn\ 



(l-V qX o(q)) 2 



(36) 



^ 2 - 



xo(q) 
i-v,xo(?) 



Using this leads to a modified result for the supercon- 
ducting coupling constant Asc- For a specific process of 
wave vector q, coupling points on the Fermi surface, we 




Vgxo(q) 

where A?„ is the unrenormalized contribution. All such 
contributions must be summed over to determine the to- 
tal Asc- The first and second factors tend to decrease 
and increase the electron-phonon coupling respectively. 
Physically, the first factor is due the screening of the 
electron-ion interaction due to conduction electrons. The 
second factor is due to the softening of particular modes 
due to the Kohn Anomaly which will in turn enhance the 
overall electron-phonon coupling. Since these renormal- 
ization factors depend on temperature through the sus- 
ceptibility xoi Tgc must be determined self-consistcntly. 

In addition to the renormalization of the Coulomb ver- 
tex, there is also the direct Coulomb repulsion between 
electrons that is taken into account through the Coulomb 
pseudopotential /ig C which is included in McMillan's 
expression**"* 



sc 



1.2 



exp 



1.04(1 + A sc ) 



where 



A*SC — 



[Asc - MscC 1 + 0.62A SC )] 



Msc 



1 



(38) 



(39) 



and use is the screened Coulomb interaction averaged 
over the Fermi surface. 

We will now estimate ^g C . Taking into account screen- 
ing within the RPA one finds 



q 1 " V qX0 (q) 



(40) 



for the screened Coulomb interaction. In Id for q w 2kp, 
V q s 0. This is due to the divergence of Xo(q) at <Z = 2k F . 
Also, one finds that for q « 0, V s (q) » l/z/(0). Using this 
RPA screened Coulomb interaction we find for our three 
band system of the (5,0) CNT 



A*sc 



A- y 



5(e kT )6(s k , T >)V T s T ,(k-k')=0.25. 



(41) 

Then, using Eq. ffity. we obtain ^g C = 0.19 for the 
Coulomb pseudopotential with the calculated values of 
the Fermi energy and Debye frequency. 

We now see how taking into account the Coulomb in- 
teraction in this manner modifies the superconducting 
transition temperature for the (5,0) CNT. The most sig- 
nificant renormalization of the total superconducting Asc 
given by Eq. (|37|l will be for the 2k F process that cou- 
ples to the out-of-plane optical mode which was previ- 
ously seen to have the overall strongest coupling. That 
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is, at temperatures where the renormalized Asc will start 
to differ from the bare Ag C , all of the renormalization 
will come from this mode. Using Eqns. 1)37(1 and 1381 we 
find that a self-consistent solution for the superconduct- 
ing transition temperature occurs at Tsc ~ 1.1 K which 
is larger than the previously calculated CDW transition 
temperature. This therefore shows that the Coulomb in- 
teractions can favor superconductivity over the CDW in- 
stability. 

For the (6,0) CNT, we see that the value of T sc with- 
out the inclusion Coulomb interaction is smaller than 
Tcdw = that we computed in the previous section 
with the inclusion of the Coulomb interaction. We there- 
fore conclude that the CDW will be dominant for the 
(6,0) CNT. 



E. Summary of Coulomb effects 

In the above, we have shown that the introduction of 
the residual Coulomb interaction will lower both the SC 
and CDW transition temperatures. We also illustrated 
the possibility that the CDW instability can be sup- 
pressed so much by Coulomb interactions that SC will 
be dominant at low temperatures. However, we stress 
the difficulty of obtaining such quantitative results. In 
principle, to obtain an accurate Coulomb interaction in 
our basis of Bloch states, one needs to use the interaction 



V kk 'k"k"> = - I d 3 rd 3 r' ^(r)^(r'). 



|-0fc''(r')-0fc"'( r ) 
(42) 



where the tpk 's are Bloch state wave functions of the CNT 
which is difficult to obtain. The Coulomb interaction V q 
we used is only a rough approximation to this more real- 
istic interaction. On the other hand, the SC and CDW 
transition temperatures have exponential dependence on 
the Coulomb interaction parameters. One also has to be 
very careful not to double-count the electron-electron in- 
teraction terms taken into account in the single-particle 
energies through the Hartree term. As shown in Ap- 
pendix^ using a method in which the charge density is 
calculated self-consistently will give more accurate values 
for the phonon frequencies calculated through the FPA. 
However, there are serious difficulties with calculating 
the bare electron-phonon vertex with such a method as 
discussed in Sec. IV CI 



VI. DISCUSSION 



A. Comparison to other carbon based materials 





KC 8 


K3C60 


(5,0) CNT 


(6,0) CNT 


(5,5) CNT 


u(Q) (eV 1 ) 


0.24 a 


0.29 6 


0.16 


0.068 


0.034 


Asc 


0.21" 


0.7 C 


0.57 


0.12 


0.031 


V (eV) 


0.875 


2.4 


3.6 


1.8 


0.92 



TABLE VI: Density of states at the Fermi energy, the su- 
perconducting coupling strength, and the attractive potential 
strength for various carbon materials. Superscripts a, b, and 
c denote Refs. [3?SI |. |4fij . and [jj respectively. 



the Fermi energy yield u(0) = 0.24 (Ref. 01) and 0.29 
(Ref. [13) states/ eV / C atom for KC 8 and K 3 C 6 n_re- 
spectively. Estimates of Asc f° r these are 0.21 (Ref. |39j) 
and 0.7 (Ref. @). In the BCS theory, Asc is expressed 
in terms of the product of the electronic density of states 
at the Fermi level and the attractive pairing potential 
strength Asc = v{Q)Vm^ Now that we have the mag- 
nitude Asc and v(Q), we can extract the magnitude of 
the pairing potential for the intercalated graphene, the 
fullercnc, and the CNTs we study. The results are sum- 
marized in Table IVT1 

The following analysis will be very similar to that of 
Benedict et al. in Ref. 9 The central idea in their analysis 
is as follows. Since curvature increases the amount of 
hybridization between a and 7r states at the Fermi energy, 
the strict selection rules for phonon scattering between 
pure 7r states in graphene will be lifted. The amount of 
(j — 7r hybridization has roughly a 1/ R dependence on the 
radius of curvature, so the matrix elements and therefore 
the attractive potential due to curvature will go as 1/R 2 . 

Neglecting presence of pentagons in fullerenes, we 
write the attractive potential for the fullerene Vbail as 
the sum of contributions from that of the graphene sheet 
Va a t and that from curvature effects I4 urve 



I, nil 



flat 



(43) 



This relation enables us to obtain the value for Vcm-vc = 
1.5 eV. Now we can write the expected attractive inter- 
action for the CNT 



V tubB (R) = Va ai + V cl 



go/2 
R 



(44) 



where Rq ~ 5 A is the radius of a fullerene and the 
factor of two comes in because there is twice as much 
(7 — 7r hybridization in a fullerene as there is in a CNT 
of radius Rq~ In Fig. 1131 we show that Eq. (|44|l . which 
was calibrated by using only quantities from intercalated 
graphene and fullerenes, is consistent with the attractive 
potentials we obtain for the (5,0), (6,0), and (5,5) CNTs. 



As a consistency check, we now compare our results 
for the attractive potential due to electron-phonon cou- 
pling to established results for other carbon based solids, 
namely the intercalated graphene KCs and the carbon 
fullerene K^Cqq. Calculations of the density of states at 



B. Beyond Mean Field Theory 

One-dimensional electron-phonon systems have several 
competing instabilities and the true ground state may be 
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1 2 3 4 5 6 7 8 

R (Angstroms) 

FIG. 13: Vtubc(-R) from Eq. 14411 calibrated with parameters 
from intercalated graphene and fullerenes (solid line) com- 
pared to the attractive potentials calculated for the represen- 
tative CNTs (filled circles). 



found only by analyzing their interplay4ia£ Hence, one 
may be concerned that we use a mean-field approach to 
analyze a Id CNT. We point out that when we calculate 
the superconducting TgQ we include the interplay of the 
CDW and SC orders. That is, the effective supercon- 
ducting coupling Asc that we obtain in Eq. 1371 includes 
softening of the 2kp phonon mode. Such an approach 
is equivalent to the two parameter RG analysis used in 
Ref.ka. The mean-field transition temperature obtained 
by our method is equivalent to the coupling constants be- 
coming of the order of unity in the RG analysis. At Tsc 
electrons start to pair, but the system has strong fluc- 
tuations in the phase of the SC order parameter. The 
most important kind of fluctuations are thermally acti- 
vated phase slips, discussed originally for superconduct- 
ing wires in Refs- EKH Phase slips lead to only a gradual 
decrease of resistivity with temperatures below Tsc- 

For an incommensurate CDW, long range order may 
not appear at finite temperature either. To understand 
the physical meaning of the mean-field transition, we can 
introduce a Landau-Ginzburg formalism.— Here we con- 
centrate on the (5,0) and (6,0) CNTs which have three 
partially filled bands with Fermi points kp and kp where 
the exact relation 2kp — kp is satisfied. We introduce 
a complex order parameter ^i(x) related to the ampli- 
tude of the lattice distortion as q(x) — e 2lklpx ^\(x) + 
e~ pX ^\{x). At low temperature the free energy is 
given by F a [*i] = J dx{a\^ 1 \ 2 + b^^ 4 + c|^| 2 ). Be- 
low the mean-field transition temperature Tcdw we have 
a < and the system develops an amplitude for the or- 
der parameter ^i. The phase of however, is still 
fluctuating, leading to short range correlations for the 
CDW order (#i(x)*£(0)) oc e" x /^ T \ Even at T = 
we can have at best a quasi long-range order for \T/i due 



to the incommensurate value of 2k p. Lattice distortions 
at 2kp can be included by introducing another complex 

field $ 2 (x) that contributes e 2lkpX <Sf 2 (x) + e~ 2lkpX ^* 2 {x) 
to the distortion amplitude. The relation 2kp — kp im- 
plies that the free energy allows coupling between \&i and 
<F 2 of the form F wa [^ u ^ 2 ] = lj dx^ff^l + <F 2 ), so 
when the amplitude of is established, it will immedi- 
ately induce the amplitude for (although none of the 
fields have a long-range order). Appearance of such am- 
plitudes should lead to a pseudogap state of the system 
below Tcdw 1 ^ The dominant contribution to electrical 
conductivity in a clean system would then come from 
the Goldstone mode of the phase of the VP's, i.e. sliding 
of CDWs (Frolich mode). Any kind of disorder (e.g. im- 
purities or crystal defects), however, gives strong pinning 
of the CDW phase and suppresses collective mode con- 
tributions to transport. Therefore, we expect insulating 
behavior of the low temperature resistivity in most exper- 
imentally relevant circumstances if CDW is the dominant 
low temperature phase. 



C. Experimental Implications 

Proximity induced—^ as well as intrinsic^ super- 
conductivity has been experimentally observed in car- 
bon nanotubes. On the other hand, the CDW state, de- 
spite being endemic to quasi Id systems has never been 
reported for carbon nanotubes. As we discuss above, 
one needs to have very small carbon nanotubes to have 
electron-phonon interaction strong enough to make cither 
the CDW or the SC instabilities appear at experimentally 
relevant temperatures. In this work we address quantita- 
tively both of these instabilities. Our main conclusion is 
that when we include Coulomb interaction between elec- 
trons, the CDW instability does no appear even for the 
ultrasmall nanotubes, whereas the superconducting Tsc 
may be in the few Kelvin range. 

In the work by Kociak et al. in Ref . |jj electronic trans- 
port through ropes of single- walled CNTs suspended be- 
tween normal metal contacts was measured. The ropes 
are composed of several hundred CNTs in parallel with 
diameters of the order 1.4 nm. It was found that be- 
low 0.5 K, the resistance abruptly drops, an effect which 
is destroyed by the application of an external magnetic 
field of order 1 T. The largest radius CNT we study is 
the (5,5) CNT, which was seen to be in the regime where 
zone-folding is applicable. For this CNT, we calculated 
-^sc = 0.031, a value far too small to support supercon- 
ductivity at this temperature even without the inclusion 
of the Coulomb interaction. This small value of Asc is 
consistent with the experimental measurements of the 
electron-phonon couplin g in CNTs of similar diameter 
by Hertel et al. in Ref. 49. It is possible that the in- 
teractions between CNTs in the rope play a tantamount 
role for superconductivity in the experiment of Ref. as 
suggested by Gonzalez in Ref. Another possibility 
is that a small number of nanotubes in the rope have 
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a small diameter. For nanotubes with a diameter of 4 
A we find superconducting Tgc in the 1 K range which 
would be consistent with these experiments. A small 
number of superconducting nanotubes could provide a 
short-circuiting in transport measurements or even in- 
duce superconductivity in other CNTs via the proximity 
effect. 

In the experimental work of Tang et al. in Ref . 0, elec- 
trical transport was measured through a zeolite matrix 
containing single- walled CNTs. In the zeolite matrix, 
the CNTs are well-separated from each other creating 
an idealized one-dimensional system. The diameters of 
the CNTs were determined to be approximately 4 A by 
measuring the radial breathing phonon mode frequency 
by Raman spectroscopy. The superconducting transition 
temperature for this system was found to be 15 K from 
transport measurements. In addition, the Meissner ef- 
fect was observed through the temperature dependence of 
the magnetic susceptibility suggesting that the large cur- 
rents observed in transport measurements are not from 
the sliding charge-density wave collective mode, but are 
indeed from superconducting correlations. 

The ultrasmall (5,0) CNT we study is the likely candi- 
date structure for the CNTs confined in the zeolite ma- 
trix in these experiments. We find for this system that 
the electron-phonon coupling is very strong. We find in 
the mean-field theory, neglecting Coulomb interactions, 
that Tcdw = 160 K and Tgc = 64 K, indicating that the 
charge-density wave instability is stronger in this approx- 
imation. However, putting in the Coulomb interaction as 
in Eq. I|28[l . the charge-density wave transition was sup- 
pressed to very low temperatures, making superconduc- 
tivity dominant with Tgc = 1 K. Discrepancy between 
our calculated Tgc and the experimentally observed 15 K 
should not be a reason for concern. The superconducting 
transition temperature in Eq. 1381 is exponentially sensi- 
tive to the strength of the Coulomb interaction, and our 
estimates of the latter are not very accurate. 



This method is elaborated in Sec. Ill Dl After the calcula- 
tion of these quantities, the effective Frolich Hamiltonian 
has been fully constructed. The remarkable agreement 
of the coefficients of the logarithmic divergences com- 
puted by using quantities from the band structure and 
the electron-phonon coupling with the frozen-phonon fre- 
quencies is a consistency check for this method. 



With the Frolich Hamiltonian, we then used the RPA 
analysis of the Peierls instability (in Sec. lIV Al ) and the 
McMillan equation (in Sec. II V Bl and Appendix |B|) to ob- 
tain the charge-density wave and superconducting tran- 
sition temperatures, the result with the higher transi- 
tion temperature being the dominant phase at low tem- 
peratures. For instance, when the CDW is dominant, 
the Fermi surface will be destroyed around Tcdw elim- 
inating superconductivity altogether. By this method, 
we provided an exhaustive analysis of three types of 
CNTs: (5,0), (6,0), and (5,5). The more conventional 
larger-radius (5,5) CNT was seen to be stable against the 
CDW and SC transitions down to very low temperatures 
(<C IK) if we only include electron-phonon interactions. 
For the ultrasmall radius (5,0) and (6,0) CNTs, however, 
the CDW was found to be the dominant phase, with tran- 
sition temperatures of 160 and 6 K respectively. For both 
of these CNTs, 2k p is incommensurate with the under- 
lying lattice. Furthermore, in contrast to larger radius 
CNTs which have dominant electron-phonon coupling to 
the in-plane phonon modes, the ultrasmall (5,0) and (6,0) 
CNTs were found to have dominant coupling to the out- 
of-plane phonon modes (see Fig. |HJ), as seen from the di- 
rect computation of the electron-phonon matrix elements 
Mkrk'r'fj,- This is further supported by the frozen-phonon 
computation of frequencies which show the most robust 
Kohn anomalies for these modes (see Fig.|6j. 



VII. SUMMARY AND CONCLUSIONS 

In this work, we have used the Frolich Hamiltonian 
written in Eq. Q to study three types of small-radius 
CNTs. For this Hamiltonian, the band structure ener- 
gies were computed by using an empirical tight-binding 
method^ to first relax the structure, and then to com- 
pute the eigenvalues of the secular tight-binding equa- 
tion. The electron-phonon interaction gurk'T 1 ^ is eval- 
uated for scattering between all Fermi points. The 
dressed phonon frequencies f2 gM are computed by us- 
ing the frozen-phonon approximation given in Eq. © 
by the displacement vectors from the dynamical matrix 
of graphene given in Ref. |25l The undressed frequencies 
17°^, which enter the Frolich Hamiltonian in Eq. JQ), are 
then extracted by using the previously computed quanti- 
ties of the band structure and the electron-phonon cou- 
pling, and the RPA analysis of the Peierls instability. 



When we include the Coulomb interaction, for the 
(5,0) CNT we find that the CDW order is suppressed 
much more strongly than superconductivity. More specif- 
ically, our analysis presented in Sec. [V] shows that the 
CDW transition is pushed down to unobservably low 
temperatures, whereas the superconducting Tgc is re- 
duced to 1 K. Hence our calculation supports the possi- 
bility of observing superconductivity in ultrasmall CNTs. 
It is quite foreseeable that a more detailed model for 
the Coulomb interaction could raise Tgc to the value 
seen experimentally, especially considering the exponen- 
tial dependence of the superconducting transition tem- 
perature on the Coulomb interaction strength. For the 
(6,0) CNT, we found that the CDW remains dominant 
when the Coulomb interactions are included due to the 
weak Coulomb interaction between electrons at the Fermi 
points, and occurs at around Tcdw =5K. 
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APPENDIX A: THE ELECTRON-PHONON 
COUPLING VERTICES 

The electron-phonon coupling matrix is given by 

dV 

M kTk , T ^ = N(i/j kr \ V — — • e qll {i)\^y T/ ). (Al) 

l 

One can see that the above expression can be evaluated 
by using the finite difference formula 

M feTfe , T> = i<iM(V«" - V )\^ k , T ,). (A2) 
u 

A method for calculating this expression with a plane- 
wave basis set was previously developed^* This section 
will be devoted to describing how to calculate M krk i T r M 
with a tight-binding method. We introduce the standard 
tight-binding notation 

\ipkr) = },A kT u\xkil) (A3) 
ii 

\Xku) = ^=Y. elk ' Rn \^i)- (A4) 

* n 

Here n runs over unit cells and i runs over basis vectors in 
the unit cell and / over orbital type. Because the kinetic 
energy operator will be the same in the perturbed and 
unperturbed Hamiltonians, we can write 

M kTk , T ^ = -(ip kT \(H^ - e F )\ipv T ,). (A5) 
u 

The reason why we keep the e F term which clearly is zero 
through orthogonality will become clear below. Expand- 
ing the wave functions in the tight-binding basis set, we 
obtain 

MkrUr>n = \Y1 A *kru{Xkii\{H (! ' 1 - e F )\xk'vv)A k , T ,in-. 

ili'V 

(A6) 

Now, we write \x q ka ) = \Xkil) + \^X k u) where the orbitals 
of \x k ii) are centered on the perturbed lattice. Inserting 
this into the above equation, we obtain 

Mkrk'r'fi = \ Y. A kril({xZ\{H w -e F )\xl^ v ) (A7) 

Hi' I' 

- {{5xZ\{H^-e F )\xk' VV ) + ^.o.))Ay Tli . v . 

In the second term we can do the substitution H qfl — > Hq 
because the effect of doing this will be second order in u 
and we are interested in an expression that is accurate to 
first order. Then, this term will be 

^2 A krti (( s Xkii\(H - e F )\xk<i'v) + h.c.) Ak'T'i'i' 

ili'V 

= Y,Al Til {5xku\(H-e F )\i> k , T ,) + h.c. = 0. (A8) 




FIG. 14: Migdal's expression for the electronic self-energy. 
The thick line denotes the dressed electronic Green's function 
and the wavy line denotes the phonon Green's functions. 



So we finally have the expression 

M kTk , T ,, = i £ A* kTU (xZ\(H^ - s F )\xl% n ,)A k , T , vl , 

ili'V 

This expression can be computed by evaluating the 
tight-binding Hamiltonian and overlap matrices for the 
distorted lattice, evaluating the coefficients A kr u and 
A k i r iiiy of the wave functions for the undistorted lattice, 
and performing the above sum. There is a slight tech- 
nical problem with the above method because k and k' 
are not the same in the tight-binding matrix. However, 
it can be shown that the correct result will be obtained 
by using (xZ\H^\x q k tn') and (xZlxZ,) for the tight- 
binding and overlap matrices (or the similar expression 
with k — > k') in the limit of a large supercell. That is, 
when the distance over which neighboring atoms inter- 
act is small compared to the length of a unit cell, this 
method becomes exact. When we apply this method, we 
checked for convergence of the coupling as a function of 
the unit cell size. 



APPENDIX B: ISOTROPIC ELIASHBERG 
EQUATIONS IN ID 

Obtaining quantitative parameters of superconductors 
described by the BCS theory like the transition tempera- 
ture and the wave vector-dependent superconducting gap 
from microscopic models has developed into a powerful 
tool for understanding experimentally realized systems as 
well as even predicting new superconductors^ Though 
excellent review articles exist^^i we will establish the 
key results of the theory below in attempt to be as self- 
contained as possible. We will also show how to incor- 
porate the electron-phonon coupling into the phonon pa- 
rameters which become important in Id due to the CDW 
instability. 

In the following to simplify notation, we will consider 
a single band system only. The central ingredient which, 
in principle, allows one to calculate the superconduct- 
ing transition temperature to high accuracy is Migdal's 
theorem^, which allows one to evaluate the electron self- 
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energy with small error as 



1 



S(fc,iw„) = -- Efc'n'p r 3 G(fc',iw„/)r3|gfcfc' M | 

x D Qfl (k-k',n-n'). (Bl) 

This expression for the self-energy is shown in Fig. 1141 In 
this equation, (3 is inverse temperature, Tj are Pauli ma- 
trices (i = gives the identity matrix while i = 1,2,3 give 
the x,y,z Pauli matrices respectively), Dq is the nonin- 
teracting phonon Green's function, and uj n = 7r(2n+l)//3 
are the fermionic Matsubara frequencies. The electronic 
Nambu-Grccn's function, a 2 x 2 matrix, is given by 
G(k,iuj n ) = (iuj n T - e k r 3 - S) _1 . 



Now, we can expand X in terms of Pauli matrices 



£ = (1 - Z)iuj n T + <jyr\. 



(B2) 



We did not include the T3 term because this just shifts 
the quasiparticle energies and similarly we neglected the 
T2 term which can be eliminated by a proper choice of 
phase for 4>. Written in terms of these parameters, the 
Green's function becomes 



G = — 



Ziu n T + E k T 3 + 4>Ti 



(B3) 



Inserting this into IB II we obtain 



x \gkk' D 0ll (k - k\n - ri). 



(B4) 



Now, we insert the identity J deS(e — e' k ) into the above 
expression to obtain 



^ 1 f , st .Z'iu) n >r +eT3-4>'ri 
£ = - / de 8(e - e k ) 



fe'n'jU 



x \gkk'^\ 2 D 0t ,(k - k', n-n'). 



(B5) 



sum over k. In the right-hand side of IB6I we then re- 
place Z(k',n') and <p(k',n') with their Fermi-surface av- 
erages Z n i and 4> n ' . This approximation is valid when 
the Fermi-surface is fairly isotropic. Now by equating 
the coefficients of the matrices To and t\ we finally arrive 
at the equations 



Z n = 1 



fnSn ^2 A(n 



(B7) 



(B8) 



where /„ = l/\2n + 1|, s„ = sgn(2n + 1), A„ = 4> n /Z n , 
and 



\{n-n') = - 



1 



^2 S{e k )S{ek')\gkk' t ,\ 2 D 0lJ ,{k~k / , n-n'). 



kk' fjt 



The Coulomb pseudopotential fig C was inserted to ac- 
count for the bare electron-electron interaction that is 
not included in our original Hamiltonian Eq. (Ip. The 
superconducting transition temperature Tsc is the tem- 
perature at which nontrivial solutions for the gap A„ be- 
gin to appear. Equations l|B7|l . JBSJ, and (|B9|) are known 
as the isotropic Eliashberg equations^ Input parameters 
have been calculated and the Eliashberg equations have 
been solved to calculate Tsc for a variety of supercon- 
ductors described by the BCS theory. We also note that 
a generalization to the case where the Fermi surface is 
anisotropic is straightforward^ 

Now, for typical three-dimensional solids the phonon 
frequencies are affected very little by the electron-phonon 
coupling. Therefore, the above formalism where we have 
used the non-interacting phonon Green's function Dq^ 
works remarkably well in 3d. This is not the case, how- 
ever, in Id where one is encountered with the CDW in- 
stability. A more accurate phonon Green's function is 
given by 



The Lorentzian term in the integrand peaks very strongly 
at e = sp — with width on the order of temperature. 
Assuming that the rest of the integrand doesn't vary as 
rapidly about e = 0, we can replace S(e-ek') with 8(ek>) 
and perform the e integral to obtain 



Z'iuj n ,T - ip'ri , 2 . , , 

==\gkW l A D 0fl {k-k ,n-nj 



k n fji V v - • 

This approximation can be seen to break down for small 
momentum (forward) scattering due to acoustic phonons. 
This case will be discussed in Appendix |0 When close 
to Tsc, <fi' wu l be small and can be neglected in the de- 
nominator of IB 61 

Now we perform the so-called isotropic approximation. 
Multiply both sides of|B6]by 8(e k )/v a (ti) where v a (Q) 
is the density of states at the Fermi level per spin and 



2n° 



(iV n ) 2 ~ (fig 



V 2 
1 1 



(BIO) 



where 0° is the undressed frequency (without electron- 
phonon coupling) and f2 gM is the dressed frequency 
(which, as seen above can have strong temperature de- 
pendence). Calculating the dressed phonon Green's func- 
tion can be challenging because one needs both fl®^ 
and Q qfl . However, we notice that when we substitute 
Eq. l|B10Jl into Eq. (|B9J) we have the fortuitous cancel- 
lation of f2° M in the numerator of D^(k,n) with that in 
the denominator of \gkk'fi\ 2 - Thus one sees that knowl- 
edge of the undressed frequencies (which are significantly 
more difficult to obtain) will not be necessary to con- 
struct the Eliashberg equations. By doing the substi- 
tution 17° 



fL„ in Eqns. i|B7l 



and ljB9|) . one 

can thereby construct the "dressed" Eliashberg equations 
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which takes into account the influence of the electron- 
phonon coupling on the phonon frequencies which is im- 
portant in Id- 
Note also that since some modes will have temperature 
dependence, the Eliashberg equations must be solved self 
consistently. That is, we must find a temperature such 
that the SC transition temperature determined from the 
Eliashberg equations is the same as the temperature used 
for the input dressed phonon frequencies. This can be 
done by iteration. Furthermore, this method allows us to 
tell which will be the dominant phase at low temperature 
of our system. If we find a self-consistent solution of the 
Eliashberg equations and Tsc > 2cdw> then supercon- 
ductivity will be the dominant correlation. Otherwise, 
the system will prefer the CDW state. 

Finally, we will write down an expression which ap- 
proximately solves the Eliashberg equations, originally 
developed by McMillan 



This integral can be evaluated to give 



Tsc = 



1.20 



exp 



1.04(1 + Asa) 



Asc - /^c(l + 0.62A SC ) 



(Bll) 



where Asc = A(0). From the above analysis, we see that 
to be self-consistent, one should use the dressed frequen- 
cies to evaluate Asc- 



APPENDIX C: INCORPORATING q^O 
SCATTERING FROM ACOUSTIC PHONONS IN 
THE ELIASHBERG EQUATIONS. 

In this Appendix, we discuss in detail the role of acous- 
tic phonons for small-radius nanotubes. Earlier theoreti- 
cal analysis of the electron-phonon interactions in Id sys- 
tems suggested that acoustic phonons can play a domi- 
nant role in stabilizing the superconducting stated*! 6 - 
We will show, however, that since the dominant coupling 
comes from optical modes, that this effect is not impor- 
tant for the CNTs we study. 

We now consider explicitly the contributions to q w 
scattering processes coupled to acoustic phonon modes 
which are not accounted for in the approximations lead- 
ing to lB6l For the electron-phonon coupling to acoustic 
modes, we take 



i\q\/L 
i + {q/qof 



(Cl) 



where go is a cut-off of order k F and L is the system 
length. We also take Sl qtl — c\q\ and Ek = v F (\k\ — k F ). 
Inserting these quantities into lBll setting Z = 1 for sim- 
plicity, we obtain for the off-diagonal element 



A (q~0) 



dq 



a*-? 2tt 

TV 

7k] 

l + (?/9o) 2 ' K-av) 2 + (ck|) 2 



2c\q\ 



(C2) 



M) _ _7 



— T- 



v F (3 *-f \ui n >\c + \uj n - w n i \v F 

qo qo 



(C3) 



\u n - uv|/c+ qo \Lu n - Lu n <\/v F + q ' 

One then sees that scattering from q w acoustic 
phonons gives an approximate contribution to Asc (when 

n = n') of Ag^T ^ = "f/(nv F c). 

Now we consider the q ~ 2k F scattering process from 
the same acoustic phonon. For this process we obtain 

n' 

I d q^^7~^- (C4) 

2-k J H u 2 n , + {v Fq y 

This integral can be evaluated to give 



A-fck 2 
f3v F 



+ (c2k F y 



-(C5) 



One then 



finds that this gives a contribution of 

^Sc^ F> = l/{^v F c) to Asc which is exactly the same 
as the q w scattering contribution. Thus one sees 
that q w scattering from acoustic phonons can be 
very important in one-dimensional electron-phonon sys- 
tems. From such a process the so-called Wentzel-Bardeen 
instability25i2£i21 can occur which has recently been stud- 
ied in the context of CNTs. 56 We also note that a similar 
analysis can be carried out for the optical phonons, and 
it is found that the q w processes are much smaller 
than the q ss 2k F process. 

With the above method, we now see how to include 
the contribution from q sa scattering into Asc- To do 
this, we simply double the contributions to Asc from 2k F 
processes which couple to acoustic phonons to include the 
5«0 contribution. In practice, we find that using this 
procedure actually changes Asc by only a small amount. 
For instance, for the (5,0) CNT, Asc only increases by less 
than 1%. This is because the dominant contributions to 
Asc are from coupling to the optical modes as discussed 
in Sec. EH 

We also point out that presence of the Wentzel- 
Bardeen singularity would significantly renormalize the 
acoustic phonon mode frequencies of the CNTs. The fact 
that the calculated phonon frequencies using the frozen- 
phonon approximation for the CNTs are quantitatively 
similar to the analogous modes of graphene as shown in 
Figs. [7| and 1111 further supports the the notion that 
the Wentzel-Bardeen instability is unimportant in these 
systems. 



APPENDIX D: LIMITATIONS OF NON 
SELF-CONSISTENT METHOD 

In this Appendix, we will discuss the limitations of us- 
ing a method in which the charge density is not evaluated 
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self-consistently. For simplicity, we will neglect the con- 
tribution from the exchange-correlation energy £xc m 
the Kohn Sham energy functional. 

First we will consider the case of the equilibrium lattice 
structure. For this, the self-consistent total energy is 
given by 




+ E cq . (Di) 

1 ion— ion \^^J 

where the charge-density is given by n(r) = J2i l^iOOl 2 ) 
Vj^jJ is the ionic potential, and E^ n - is the ion-ion in- 
teraction. In the above and in what follows, the i summa- 
tion is carried out only over occupied electronic states. 
Applying the variational principle to Eq. ID1I gives the 
equation for the wave functions and therefore the 
charge-density n(r) 

H°>]|^) =£i|Vi> (D2) 

where 

n ° q[n] = £ + v ™ ir) + 1 dV F^r (D3) 

In solving this equation, the charge density n(r) enter- 
ing H cq [n] must be determined self-consistently to agree 
with the eigenfunctions tpi. Using this, the self-consistent 
totally energy for the equilibrium lattice is determined to 
be 

E t% = £>iiw eq [*> + pcq w (° 4 ) 

i 

where 

F c «[n] = ~\J dW 7 ^! + EZ_ ioa . (D5) 

to be essentially the same as for non- interacting atoms. 
In the tight-binding limit we expect the equilibrium 
electron density to be essentially the same as for non- 
interacting atoms. If we denote the latter as rc.o(r), we 
can replace n(r) by no(r) in Eg^ and expect the resulting 
non-self-consistent total energy -E^sc t° ^ e Quite close to 
the self-consistent total energy for the equilibrium lattice 
structure: 

EffscKEg.. (D6) 

This approach is the basis for using an effective tight- 
binding model to calculate band structures. 

Such a method, however, breaks down when we con- 
sider a lattice perturbed by a phonon. In the presence 
of a lattice distortion, the ionic potential changes to 
VJon* = V^ q + SVi on which, in turn, makes the charge- 
density non-uniform n = uq + Sn. The energy of the 
distorted structure is then 

*& = B*l + + \ J hW 

+ Eg?.*,. (D7) 



Now replacing n with hq + Sn, this can be written as 

^SC = ^(^|W dist [n ]|^)+F dist [n ] (D8) 

i 

+ \ J d\d?r< (y 1 ^ Sn(r)Sn(r') 

where H dist and F dist are given by H cq and F cq de- 
fined above with V|° q and E°^ n _ ion replaced by VJ d „ st and 

-^ion-ion- The nrst tw0 terms 011 the ri S ht of E q- EH 
can be seen to be the total energy of the distorted struc- 
ture computed with the non-self-consistent method. We 
therefore obtain 

El¥ = + \f d 3 rd 3 r' (j^^TJ ) Sn(r)Sn(r'). 

(D9) 

Subtracting Eq. ID6I from this then gives 

AE SC = A£ N sc + \ J d 3 rd 3 r' \r^i\) 5n{v)5n{r') 

(D10) 

where A-E S C,NSC = #sc!nsc _ S SC,nsc- Rewriting the 
second term in momentum space gives 

AE SC = AE NSC + l -J J^Lv(q)\5 n(l \ 2 . (Dll) 

which then makes it clear that AEsc > Ai?NSC- So 
we see that using the a non-self-consistent method to 
calculate phonon frequencies by the frozen-phonon ap- 
proximation will underestimate the phonon frequencies. 
More specifically, in a non-self-consistent method, the 
Hartree term displayed in Eq. 1D11|) is not accounted for. 
This should be particularly important in the vicinity of 
a CDW instability, where there will be a larger response 
of the charge distribution to a lattice distortion. 

APPENDIX E: DERIVATION OF EQ. (|3ljl 

In this Appendix, we will derive Eq. I|31[) by evalu- 
ating the integral appearing in Eq. I|3U|) . To estimate 
this Coulomb interaction integral, we will take the tight- 
binding wave function of graphene 

(El) 

Now k is a two dimensional vector in reciprocal space 
of the graphene lattice and 7 = ±1 corresponds to the 
conduction and valence bands. Orbitals centered on the 
first and second carbon atoms respectively in the nth 
unit cell are given by </>„i(r) and n2 (r) respectively, and 
/(k) is given by f(k) = 1 + e ~ 4kai + e ~* k ' a2 where ai 
and a 2 are the lattice vectors of graphene. For metallic 
large radius CNTs, the Fermi points correspond to K = 
■|(bi — b 2 ) and K' = |(bi — b 2 ) where bi and b 2 are 
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the reciprocal lattice vectors corresponding to ai and a.2 . 
For these points, we have /(K) = /(K 7 ) = 0. However, 
for the smaller radius CNTs we study, as indicated by 
the failure of the zone-folding method, the Fermi points 
are shifted away from K and K'. We denote the Fermi 
points of the inner band r a of the (5,0) CNT by k Ta + = 
K + k x x — k y y and k Ta _ = K — k x ~k — k y y and for the 
other inner band Tb by k Tb+ = K' + k x 5t+k y y and k Tb _ = 
K' — k. x x + k y y where the x-direction is still along the 
CNT axis and the y-direction is along the perimeter. 

For backward scattering, we take q rs 2kp , k s» 
— fcp, rs kp. Keeping only products of Carbon orbitals 
centered on the same atom, we obtain 



'-|0»i(r)| 2 



N^" 2 

n 

/*(k Ta+ )/(k T „_) 



i/*(k Ta+ )/(k To _: 

+ |<Mr)| 2 - 



Now we make use of the slow variation of e lqx compared 
to the localized orbitals to write 



-iqx}_ \ 1 

N ^ 2 

n 

r(k Ta+ )/(k ra _) 
i/*(k Ta+ )/(k To _)i 



(E3) 



I0m(r)| 2 



where t = -4ox is the basis vector for the second Carbon 

v3 

atom in the primitive unit cell. Finally, in evaluating the 
integral in Eq. Q30[) it is sufficient to replace the functions 
Jr J2 n I0ni,2(r)| 2 which vary more rapidly than V(r) by 



their average values. That is, we substitute 

/*(k ra+ )/(k Ta _) 



V'fc +gr (r)V'fer(r) 
1 



(E4) 



■iqx 



2ttRL 



JqSit 



|/*(k Ta + )/(k Ta _ 



Using the same approximations for the factor 
4 } k'-q{ Y ')4 } k'{ Y '): w e obtain for the Coulomb inter- 
action 



1 

4 
1 

1? 



/*(k ro+ )/(k Ta _) 



Jqx-t 



|/*(k Ta+ )/(k To _)| 



(E5) 



dxdx'e- tq{x - 



2ttR 



dy 

2vri? 



2ttR 



dy' 
2nR 



V(r - r'). (E6) 



We will now evaluate the prefactor in this equation 
for the inner band of the (5,0) and (6,0) CNTs. Using 
the calculated Fermi points along with the zone-folding 



. 0.11 2tt 

(5,0) CNT. From this we obtain 



method, we obtain k x = ±^4i— and k v = rk— for the 

' x ^/3 a y 15 a 



r(k Ta+ )/(k Ta _) 



JqSt-t 



l/*(k Ta+ )/(k Ta _; 



= 0.59 



(E7) 



For the (6,0) CNT the Fermi points are k x = ± — — 
and k y = 0. This gives 



0.076 2w 



/*(k Ta+ )/(k Ta _) + el9 -. t 



l/*(k Ta+ )/(k Ta _) 



0.0016 



(E8) 



which is smaller due to the different symmetry of the 
wavefunctions at the Fermi points. These are the values 
of the prefactor 7 appearing in Eq. i|31|) . 
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